On the Nature of Infrared Singularities in d < 2 Disordered Interacting Systems. 
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We address the problem of infrared singularities in the perturbation theory for disordered interacting 
systems in d < 2. We show that a typical, sufficiently large interacting system exhibits a linear 
instability in the spin triplet channel. In the density-density channel, although stability is preserved, 
a large number of soft modes is accumulated. These phenomena are responsible for the instability 
of the weak-interacting fixed point. Although generic, the unstable direction and soft modes are 
highly sample specific and can not be effectively captured by conventional techniques based on 
an averaging procedure. Rather, the instability is determined by the largest eigenvalues of the 
polarization operator. We propose to employ the optimal fluctuation method for evaluating the 
probability of such events. 

PACS numbers: 71.10.Hf, 71.30. +h, 73.23.Ps 
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The physics of two-dimensional {d = 2) disordered 
electronic systems has received a new impetus since 
the accumulation of experimental evidence of a metal- 
insulator transition in zero magnetic field [Q. Other re- 
markable experimental developments are associated with 
d = 2 finite size structures - quantum dots . The data 
indicate the crucial role of e-e interactions in these disor- 
dered (chaotic) systems. We also point out that in both 
cases spin degrees of freedom seem to play a major role. 

On the theoretical side it was known since the work 
of Altshuler and Aronov |^ that various physical observ- 
ables acquire infrared singular corrections in the presence 
of e-e interactions. Using a finite system size, L, as a cut- 
off one may write the correction to an observable, x 
(which may be specific heat, conductivity, spin polariza- 
tion etc.) as 



Sx L^-'^ 
— (X 



(1) 



where Pd is a d-dimensional density of states at a Fermi 
level and D is a diffusion constant. Hereafter, L^^^ 
should be understood as In L /I with I being a microscopic 
length (mean free path ||^) for the d = 2 case. Once the 
system size becomes comparable with the characteristic 
length = {ydDY'^^'"^^ (6 = le"^"), the validity of the 
first order perturbation theory breaks down, and expres- 
sion (|l|) can not be applied. Despite the fact that has 
an interpretation as the localization length, Eq. (|l|) has 
nothing to do with the actual Anderson localization. In 
fact, in the presence of a small magnetic field the physical 
localization length may be much larger than S^d- Thus, 
the behavior of an observable x at L > may be studied 
quite independently of the onset of the Anderson local- 
ization. 

Finkel'stein |^,^ has developed a renormalization 
group approach to the problem. His results show that 
a weak-disorder and weak-interaction fixed point is un- 
stable at d < 2. As a result, the large scale behavior of 
the system is governed by some other fixed point. He has 



also noticed that the system acquires a tendency towards 
a partial spin polarization and pointed out that the spin 
susceptibility diverges independently of all (and before) 
other quantities 0. In spite of considerable further ef- 
forts in this direction [p|-[Tot, the physical nature of the 
new fixed point remains to a large extent unclear. The 
ultimate fate of the corrections, Eq. (|l|), in the region 
L > is also unknown. 

The purpose of this letter is to shed some light on the 
above questions. We concentrate primarily on the spin 
triplet interaction channel, restricting ourselves to few re- 
marks concerning spin singlet (density-density) channel. 
We argue that for each particular realization of disorder 
and L > ^d the system exhibits a linear response in- 
stability in the spin triplet channel. On the mean-field 
level this instability leads to the spontaneous breaking 
of the spin-rotation symmetry and partial spin polar- 
ization. The questions concerning stability of true long- 
range order with respect to fluctuations in the thermo- 
dynamic limit are well beyond the scope of this letter. 
For a finite size system at sufficiently low temperatures, 
however, the fluctuations will not destroy spontaneous 
magnetization. The existence of an unstable direction 
(in the functional space of possible spin polarizations) is 
a generic property of any disordered system. Its concrete 
shape, however, depends on the polarization operator for 
a given disorder realization and is extremely sample de- 
pendent. Methods which employ an ensemble averaging 
procedure at an early stage of the calculations are bound 
to loose information about this instability since they use 
the disorder-averaged polarization operator. Instead of 
being linear, the effect appears to be encoded into higher 
non-linear corrections. The stability of the paramagnetic 
state for a given system has very little to do with the av- 
eraged susceptibility. The importance of fluctuations of 
various susceptibilities in disordered metals was stressed 
by several authors ||ll|,|l^ (see also |jl^ for applications to 
superconductors). In fact, the stability of the paramag- 
netic state is determined by atypically large eigenvalues 
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of the polarization operator |T^. Thus, methods treat- 
ing the optimal fluctuations seem to be more adequate. 
These methods were initially developed for the problem 
of "Lifshitz tails" by Zittartz and Langer and used re- 
cently in the problem of pre-localized states . Below 
we adopt them for the treatment of interacting electrons. 

Consider a gas of interacting electrons, moving in a 
random potential. In terms of the slow degrees of free- 
dom the interaction Hamiltonian takes the form 



Vs{r — r')n{r)n{r') — Vts{r)s{r) 



(2) 



where the charge, n{r) = ^ ip'jj{r)'ipa{f')^ and spin, s(r) — 

a 

^ 'i/'J.(r)cr(j(j''0(T'(?'), densities are slowly varying on the 

crcr' 

scale \f\ Vt = VtS{r — r'). We write the partition func- 
tion as an imaginary time functional integral [ p^ and 
perform the Hubbard-Stratonovich transformation of the 
singlet and triplet interaction terms by means of fields 
and H correspondingly. Having in mind developing a 
Landau-Ginzburg mean-field theory at finite tempera- 
ture T, we restrict ourselves to the zero Matsubara fre- 
quency only, H(r) = H(r, = 0). The subsequent 
integration over the fermionic degrees of freedom results 
in det[l + G(i$ + Her)], which may be now expanded in 
powers of $ and H. In this way one obtains 



(3) 



where fH^ designate non-linear terms. The central 
quantity of interest is the static polarization operator 
(PO), defined as 



(4) 



where Gi^{r,r') is a sample specific Green function and 
£„ = 7rr(2n+ 1). The Gaussian part of the functional in- 
tegral, Eq. (^, becomes unstable if the positively defined 
Hermitian operator II(r, r') has at least one eigenvalue 
larger than Vj~^. We are faced, thus, with the spectral 
problem for the PO: 



dVn(r,r')$„(r') - A„$„(r): 



(5) 



J |$|^(i''r — 1. Specifically, we are interested in the large 
A tail of the PO spectral density since the spontaneous 
symmetry breaking will be determined by the largest 
eigenvalue. The spectral density is defined as 



iV(A)^(^5(A-A„) 



(6) 



where the angular brackets denote disorder averaging. 
Employing the optimal fluctuation we shall demonstrate 
that 



iV (a > Atypical) occxpj-/fA 



VdD 



(7) 



with f{x) being a certain universal function (~ a;^ ~ in a 
simplest approximation). Eq. demonstrates that the 
condition for the appearance of a large (A > V^"^) eigen- 
value is L^^'^ > VdD/Vt which coincides precisely with 
the breakdown of the perturbation theory, Eq. (|l|). It is 
clear now from Eq. (|^) that once such an eigenvalue is 
formed, the system develops a non-trivial saddle point 
with H(r) 7^ 0, which describes a spin polarized state. 
A new expansion arround this sample-specific minimum 
should be developed. 

The (replicated) ensemble averaged version of the ef- 
fective action in Eq. (||) is 



V,-' - (ft) 



H 



;f)--((ft2» e 



(8) 



The average PO, (n(r, r')) = 6{r — r'), is an operator 
with all eigenvalues equal to unity and thus with no tails 
in spectral density . This fact ensures the stability of 
the Gaussian integral for Vj~^ > 1, which is the Stoner 
criterion for the ferromagnetic instability [p^ . The in- 
formation about anomalously large eigenvalues of 11 is 
now hidden in the higher order non-linear terms. Thus 
the basic fact about the presence of a linear instability 
in the theory appears to be obscured by the averaging 
procedure. 

Let us make a few remarks concerning the singlet in- 
teraction channel. The corresponding functional integral 
over a scalar field $ has a Gaussian part of the form 



^ / / d'^rdV$(r) 



V-^{r-r')+U{r,r') $(r') . (9) 



Note the plus sign in front of PO! The kernel is now 
strictly positively defined and can not exhibit an in- 
stability. It may have, however, very soft (and long- 
wavelength) modes associated with a small A (A > 0) tail 
of the spectral density A^(A) These atypically soft 

modes may lead to singularities in the singlet channel. 
One should be careful with these arguments, however, 
since there is no reason to neglect non-zero Matsubara 
frequencies in this situation. For a;,„ > DQ^ even the av- 
erage PO has small eigenvalues. The fact that the diver- 
gence of singlet and triplet quantities is associated with 
the opposite tails of A^(A) may be related to Finkel'stein's 
results quoted above. 

We turn now to the evaluation of the tails of the PO 
spectral density, N{X). To this end we define the corre- 
sponding Green function operator as 
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and A'^(A) = Tr~^QTT{Qx). To perform disorder averag- 
ing we need to know the statistics of the random operator 
n(r, r'). A straightforward diagrammatic calculation of 
its second cumulant yields at T = 



^((fl(r,r')n(.',.))) 



(11) 



edeDe{r - s')D^{s' - r')D^{r' - s)De{s - r) . 



Here D^{r) is a propagator of the classical diffusion op- 
erator 



Fd 



DQ^ + \e\ D \L 



(12) 



where = y^D/e and Fd{x) is a dimensionless func- 
tion with the assimptotic behavior -Fd(O) = const and 
Fd{x oo) ^ e^^. Performing the energy integration in 
Eq. (O) one obtains 



(13) 



We can perform now the averaging of the Green function, 
Qx, assuming the Gaussian distribution for 11 . This is cer- 
tainly not an exact procedure because the higher order 
cumulants of 11 are not negligible in the most interesting 
parameter region L > £^d- We shall follow, however, this 
idea since it allows us to illustrate the method and pro- 
vides a basis for further generalizations. We thus obtain 



4n. 



(14) 



where integration over repeated coordinates is implicit. 
The average value of the PO is omitted in this expres- 
sion, since it leads to a trivial redefinition A — > A — 1 
only. We argue now, following Zittartz and Langer [H , 
that large A tails are determined by the saddle point of 
this functional integral. Variation over 11^^' leads to the 
equation for the optimal realization of the PO 



n(r,r')=QI 



■ gx{R,s)gx{s',R') 
Qx{R,R') 



(15) 



In the close vicinity of some atypically large eigen- 
value, Ao, the Green function may be well approximated 
by the single term in the sum, Eq. (p^, Qx{fTr') = 
<^o{r)%{r')/{X - Ao), leading to 



n(r,r') = 



1 



A- Ao 



(16) 



Again the d-dimensional integration is assumed over re- 
peated coordinates. Finally, substitution of Eq. (16) into 
Eq. (^) results in the nonlinear self-consistency equation 
for the typical eigenfunction corresponding to a large 
eigenvalue, Aq, 



Ao$o(?') = 



A - A, 



■QZ,^,{t'W,{s)^,{s'). (17) 



To avoid the solution of this non-linear integral equa- 
tion one may use purely dimensional arguments. To this 
end let us write this equation in terms of dimensionless 
coordinates x — r/L and dimensionless eigenfunctions 
(t)Q{x) = L'*/2$o(r-). We also employ Eq. (|l|) for the 
cumulant Q: 



00 (x) 



L2 



(18) 



VdD^J\o{\ - Ao) 
{\x - y'\ + \y' - x'\ + \x' -y\ + \y- a;|)4('i-i) 



To ensure the existence of a normalizable 

J |0op(i''a; = 1 j solution of this equation one must 
/ 
require that the expression in the square brackets on its 

r.h.s. is of order unity. We conclude that 



(19) 



The last step is to find the statistical weight of the 
optimal realization n, which has Ao as an eigenvalue. 
Its statistical weight is obviously given by IIQ^^n — 
Ao/(A — Ao), where we employed Eqs. (p^, (17) and nor- 
malization condition for <I>o. The last expression together 
with Eq. (n9|) implies that 



iV(Ao) oc exp { -Cd Ao 



>l2-d 



(20) 



where Cd is a numerical constant of order unity. The 
fact that Gaussian tails are obtained may be traced back 
to the assumption about the Gaussian distribution of 
the PO and should not be taken too seriously in the 
regime L > S^d- However, the dependence of the spec- 
tral density on the parameter XvdD / L?~^^ follows from 
pure dimensional analysis and may be justified for more 
realistic assumptions concerning PO statistics. In fact, 
the Gaussian spectral density, Eq. (pO|), is the most pes- 
simistic estimate for the tails. Indeed, all the effects we 
have neglected (higher order cumulants of H, and onset 
of Anderson localization) should increase fluctuations of 
the PO, making the tails decrease slower than Gaussian. 
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It remains a challenging problem to calculate the spec- 
tral density of the PO under more realistic assumptions. 
The "pessimistic" result, Eq. (pO|), is already sufficient to 
demonstrate our main point: one may find an arbitrarily 
large eigenvalue of the PO if the system size is taken to 
be large enough. This leads to a ferromagnetic instability 
in the triplet channel for arbitrarily small interaction, Vt- 

In the most interesting d = 2 case the argument in 
Eq. ( pO| ) has the form Xog{hiL/l)^^, where g = V2D > 1 
is a dimensionless conductance of a 2d metal. To find an 
eigenvalue of the PO Aq > Vj~^ > 1 one should typically 
consider a sample with a size, L > leyi'p{g/Vt\. This 
still may be much smaller than the localization length 
in the unitary ensemble, which is of the order /expj^'^}. 
This is to say that, contrary to a naive interpretation of 
Eq. (po|), one may encounter the ferromagnetic instability 
well inside the metallic side of the Anderson transition. 

We want to point out that in low dimensional [d < 2) 
systems the spontaneous symmetry breaking occures into 
an extended state, rather than into a localized one (as 
it happens in 3d systems [^). This is clear from the 
extended nature the solution of Eq. (|l8|). In this case 
the role of nonlinear terms ( FH'* in Eq. (||) ) reduces 
to merely determining the amplitude of the spontaneous 
globally coherent magnetization. 

The aim of this letter is to elucidate the source of some 
problems in the theory of d < 2 disordered interacting 
electrons. The traditional perturbation theory and RG 
treatment lead to infrared divergence and instability of 
the weak-coupling fixed point correspondingly. They do 
not provide a simple physical reason for the singularities, 
nor do they explain the nature of the strong coupling 
fixed point. We argue that at least some of the difficul- 
ties can be traced to (i) the ferromagnetic instability of 
a given system in the triplet channel, and (ii) to the ac- 
cumulation of many soft modes in the singlet channel. 
The latter necessitates taking into account non-linear 
terms in the action. Both the unstable directions and 
the soft modes, although generic, are sample specific and 
do not survive (traditional) ensemble averaging proce- 
dure. It would be desirable to construct a theory, which 
first adjusts the integration directions (bosonic fields) to 
the concrete sample configuration and only then performs 
the averaging. 
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